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Abstract 

This paper presents a one-dimensional steady-state mathematical thermal power model of the ASRG. 
It aims to provide a guideline of understanding how the ASRG works and what can change its 
performance. The thermal dynamics and energy balance of the generator is explained using the thermal 
circuit of the ASRG. The Stirling convertor performance map is used to represent the convertor. How the 
convertor performance map is coupled in the thermal circuit is explained. The ASRG performance 
characteristics under i) different sink temperatures and ii) over the years of mission (YOM) are predicted 
using the one-dimensional model. Two Stirling converter control strategies, i) fixing the hot-end of 
temperature of the convertor by adjusting piston amplitude and ii) fixing the piston amplitude, were tested 
in the model. Numerical results show that the first control strategy can result in a higher system efficiency 
than the second control strategy when the ambient gets warmer or the general-purpose heat source 
(GPHS) fuel load decays over the YOM. The ASRG performance data presented in this paper doesn't 
pertain to the ASRG flight unit. Some data of the ASRG engineering unit (EU) and flight unit that are 
available in public domain are used in this paper for the purpose of numerical studies. 

1.0 Introduction 

The advanced Stirling radioisotope generator (ASRG) (Refs. 1 to 3) is being developed for 
multimission applications to provide a high-efficiency power source alternative to radioisotope 
thermoelectric generators (RTGs). The ASRG efficiency could reach 28 to 32 percent, which results in 
reducing the required amount of radioisotope by roughly a factor of 4 compared to RTGs. Thus, because 
of the limited supply of Pu-238, utilization of the ASRG can extend radioisotope power available for 
future space science missions, such as deep-space missions, large planetary surface rovers, and systems in 
support of human exploration activities. 

An overview of the ASRG is shown in Figure 1. The ASRG consists of two advanced Stirling 
convertors (ASCs) enclosed in the housing; each has a general purpose heat source (GPHS) attached at 
the hot end to provide the heat. A gas management valve (GMV) and pressure relief device (PRD) are 
located at the top of the housing. The housing with attached fins radiates the heat to the environment. The 
GMV is used to maintain a near-atmosphere pressure of inert gas inside the housing during ground 
operations. This gas is permanently vented to vacuum by the PRD when the ambient becomes vacuum. 
The ASC control unit (ACU) is separate from the ASRG housing. It converts the AC signals from both 
ASC to 28 to 34 VDC for a typical spacecraft electrical bus. The controller is used to maintain 
synchronized displacer/piston movement of the two directionally opposed Stirling convertors to minimize 
induced disturbance to the spacecraft and its precision instrumentation. 

An ASC consists of a free-piston Stirling converter and an integral linear alternator that converts the 
piston reciprocating motion to electrical power output. The Stirling engine is a heat engine that operates 
by cyclic compression and expansion of a working fluid at different temperature levels such that there is 
net conversion of heat energy to mechanical work. The GPHS provides the heat to the heater head of the 
ASC. Helium is used as the working fluid that is hermetically contained within the converter pressure 
vessel and is in the various working stages of the ASC. The displacer shuttles helium between the 
expansion space, where heat is received, and the compression space, where the heat is rejected at the cold 
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Figure 1 . — Overview of the ASRG components 


side adapter flange (CSAF) attached to the housing. The changes in pressures and volumes of the 
convertor working spaces drives a power piston that reciprocates to produce AC electrical power via a 
permanent magnet linear alternator. 

The ASRG DC power output depends on the following: 

1) GPHS fuel load at the year of mission (YOM) 

2) Thermal insulation loss from GPHS to the environment 

3) ASC efficiency (i.e., Stirling engine efficiency x alternator efficiency, the ratio of the AC power 
output and heat going into the engine) 

4) Cable/connector power loss 

5) ACU power loss 

In this paper, a one-dimensional steady-state thermal power model is built to predict the ASRG AC 
power output. The ASC is modeled using the Sage performance map (Ref. 4) that provides the ASC 
efficiency. With given GPHS fuel load and defined ambient conditions, the AC power output is computed 
using the one-dimensional model. Assuming the cable/connector power loss is 0.8 percent of the AC 
power and the ACU efficiency is 87 percent, the ASRG DC power can be computed and will be presented 
in the paper. 

In the following, the ASRG thermodynamics and energy balance will be explained first, then ASRG 
performance under varying ambient sink temperature (for different missions) and over the year of mission 
(for GPHS decay) are shown using two ASC control strategies: (i) fixing the hot-end temperature of the 
ASC and (ii) fixing the piston amplitude; which is followed by the summary and conclusions. 

2.0 Thermal Dynamics and Energy Balance of the ASRG 

2.1 The Fundamentals of the Idealized Stirling Cycle 

The idealized Stirling cycle is a Carnot cycle (see Figure 2) with four thermodynamic processes 
acting on the working fluid: 

1) Isothermal expansion at the hot-end temperature of the ASC, T h 

2) Constant-volume heat removal 

3) Isothermal compression at the cold-end temperature of the ASC, T c 

4) Constant-volume heat addition. 
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Figure 2. — P-V diagram of the Carnot cycle 


The heat going into the engine, Q in , is the area of ABFE. The heat rejected from the engine, Q reJ , is the 
area of CDEF. The mechanical work, W (= Q„, - Q rej ), is the area of ABCD. The net power output, P= W 
f with /being the frequency of the Stirling engine. It was shown in (Ref. 5) (Eq. (14)) with K= 0, T h = 7/ 
Qin = g 3 - 4 , Equation (10) with K = 0, / = 7j ,Q reJ = - Q u2 ), we can have 

Qm=mR]n Kll±^L Th (1) 

V D 

Q =mRhX°±^T c (2) 

V D 

and 

p = (Qin Qrej )/ = finR In (T h -T c ) (3) 

V D 


where V D is displacer swept volume, V P is the piston swept volume, m is the mass of the working fluid, 
and R is the gas constant. From Equations (1) and (2), then we can have the resistances at the hot and cold 
ends, respectively, of an idealized Stirling engine defined as 


ni _ Th 
h ~ 


0 s T1 1 Tn +Vp 

finR In — - 




1 


6 - C Q« 


(4) 

(5) 


However, in practice, it is difficult to determine V D and V P . Some uncertainties of the Stirling engine 
performance such as regenerator effectiveness, dead volumes, and other losses have to be considered. The 
shaded area inside ABCD in Figure 2 represents the actual work done by the converter. West (Ref. 6) gave 
an engineering model to compute the power output of a Stirling engine as follows: 

P=w n^r?-(T h -T c ) (6) 

1 h + 1 c 


where p m is effective mean pressure of the Stirling engine and w„ is West number for the Stirling engine, 
in the range of 0.25 to 0.35. 

Following Equations (4) and (5), we can define the resistance based on the West equation as follows: 
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JVpPm 

T h + T c 


It can be seen that the resistance at both the hot and cold ends are the functions of the Stirling engine 
frequency, piston swept volume (piston amplitude), effective mean pressure, T h , and T c . When the piston 
amplitude decreases, the resistance will increase. The resistances at the hot and cold ends are the same 
based on equations for an idealized Stirling engine or West equation. However, in reality, the resistance at 
the hot and cold ends of the engine could be different due to the losses in the engine are counted in 
different ways at the hot and cold ends, which is shown in the Sage performance map (Ref. 4). 


2.2 Thermal Circuit of the ASRG 

The thermal circuit for half of the ASRG is shown in Figure 3, where T s is the temperature of the GPHS, 
T h is the temperature at the hot end of the ASC, T c is the temperature at the cold end of the ASC, T w is the 
housing wall temperature and T a is the equivalent ambient sink temperature. R\ is the contact and 
conduction resistances between GPHS and housing wall; R 2 is the contact and conduction resistances 
between GPHS and heat collector of the ASC; /G, is the contact and conduction resistances between the ASC 
hot end and housing wall; R 4 is the conduction resistance between the cold end of the ASC and housing 
wall; R\ is the convection and/or radiation resistances between housing wall and the ambient environment; 
R 5 is the conduction resistance between the hot and cold ends of the ASC; R b j is the resistance of the ASC 
at the hot end. The convertor acts like a thermal resistance in the thermal perspective. 

With Q tota i and T a being defined and imposed as boundary conditions, the total heat Q, ota i provided 
from the GPHS will split into Qi oss i(loss through the thermal insulation around the GPHS) and Q in to the 
ASC. The heat going into the ASC further splits into three paths: Qi 0S s 2 (loss through conduction), Qi oss j 
(loss through insulation), and Q wor k ( heating up the gas to do the work). (Q rej is the heat rejected from the 
ASC through the CSAF attached to the housing. All the waste heat will be rejected to the ambient through 
radiation and/or convection. How the heat splits depends on the thermal resistances in the circuit, such as 
R\, R 2 , R 2 , R 4 , /G, R' 4 and R 6 h (see Figure 3). Re j, and Q rej are functions of T h , T c , Q in , Q rej , and A p (piston 
amplitude). R\, R 2 , /G, R 4 , /G, R' 4 can be defined based on the material properties and the condition of 
contact surfaces. R 6 j, and Q rej will be defined based on the ASC performance map (Ref. 4). T h is a 
boundary node used to define R 6 h between T h and T b Q rej is imposed as a heat load at a boundary node 
that has temperautre of T c An iterative procedure has to be used to get the correct R () h and Q rej that satisfy 
both the energy balance of the ASRG thermal circuit and the ASC performance map. The three steps in 
the iterative procedure are listed as follows; (i) give initial values of R b j, and Q rej ; (ii) run the thermal 
model to compute temperatures and heat flows, such as T h , T c , and Q„, : (iii) use T h , T c , and Q m to update 
the R b /, and Q re j based on the ASC performance map; The three steps are repeated until the differences of 
Re j and Q rej between two iterations are small enough (10 4 ) to conclude that the solution is converged. 

Note that housing temperature can be quite different on the housing and the convective heat transfer 
could vary significantly around housing when it exists. Here, the variation of the housing temperature is 
not considered. It is assumed that only radiation exists in the following ASRG performance prediction. 

In the following, the relationship between Q in , T h and R b j (i.e., A p ), T a at the hot end will be derived 
first, which is followed by the relationship between Q rej , T c and R 6 c (i.e., A p ), T„ at the cold end. 
Assuming Q in = Q wori k to simplify the derivation, we have 


Qtotai ~ Qioss _ 1 + Qin 

(8) 

T -T 

o — s 

x-'loss 1 n 

K \ 

(9) 
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( 10 ) 

(ii) 


Equations (10) and (11) give 


T s =Q m R 2 + T h (12) 

T h =Q in R 6 H (13) 


Substituting Equations (12) and (13) into Equation (9), we will get 

( t ? 2 +R 6 i, }Qi n ~T W 


Qloss 1 




By substituting Equation (14) into Equation (8), we will get 

T , 


Qtotal + 


Qi, 


R 


1 + 


Ro + R, 
hi 


6 h 


Then substitute Equation (15) into Equation (13), and we get 

T , 


Qtotal E 


R> 


Th 1 + R 2 /R l 1 


R 


6 h 


R< 


(14) 


(15) 


(16) 


It can be seen that if T w increases (i.e., T a increases), T h and (Q, n will increase and Qi oss \ will decrease. If 
Re h is increases, T h will increase, Q in will drop, thus Qi oss j will increase. 

At the cold end of the ASC, we have 
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(17) 


Qrej A £ V( T W ~T a 4 ) + hA 2 (T w -T a ) = 

R\ 


-re] 


T c ~ T w 

Ra 


(18) 


Q« = 



(19) 


where A\ is the housing area, A 2 is the convective heat transfer area, s is the surface emissivity, cr is the 
Stefan-Boltzmann constant, and h is the convection heat transfer coefficient. Equations (17) and (18) give 


T w~ T a + R 4Qrej 

T W =T C - tUQrcj 

By combining Equations (19), (20), and (21), we get 


T = 

C 


^ R 4 + R' 4 

Ra 


l 6 c 


T =- 


1- 


*4 
R 6 c~ R 4 


( 20 ) 

( 21 ) 


( 22 ) 


(23) 


It can been seen that when T„ increases, T c and Q ny will increase. When R b c increases, T c and Q rej will 
decrease. 

In summary, 

1 ) If the piston amplitude decreases, 
o R 6 j and R 6 _ c increase, then 

o T h increases; Equation (16) 
o Q ln drops; Equation (15) 
o T c drops; Equation (22) 
o Q rej drops; Equation (19) 
o the net power output ( P ) drops. 

2) If the sink temperature increases, 
o T h increases; Equation (16) 

o Q in increases; Equation (15) 
o T c increases; Equation (22) 
o Qrej increases; Equation ( 1 9) 
o the net power output (P) drops. 


3.0 One-Dimensional Thermal Power Model Results 

The one-dimensional thermal model described above is built in Matlab and used to predict the ASRG 
performance using two ASC control strategies. In the following, the ASRG performance characteristics 
under Y OM is described first, which is followed by the corresponding performance under different 
ambient sink temperature. 
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3.1 ASRG Performance Characteristic Under YOM 


The ASRG power output is studied for the YOM since the GPHS fuel decays with time. A 17-year 
time frame that includes up to 3 -year storage and 14-year on mission is considered, and the fuel decay 
rate is computed as 


YOM/ 

YOM GPHS fuel load = BOM fuel load x 0.5 /s7 74 

The begining of mission (BOM) fuel load is 244, 250, and 256 W for minimum, nominal, and maximum 
cases, respectively. As shown in Figure 4, for BOM fuel load = 244 W, each GPHS fuel load drops 1.8 W 
per year and has a 30.6-W total drop at the end of mission (EOM). For BOM fuel load = 256 W, each 
GPHS has a 32.0-W total drop at the EOM. The fuel almost drops linearly with time. 

With the fuel decaying over the Y OM, the ASRG power output in deep space (no sun, T a = 4 K) is 
shown in Figure 5(a) using the first ASC control strategy(fixed T h ). It can be seen that ASRG DC power 
output drops 1 9 W within the 1 7 -year period; that is, a 1 . 1 2-We drop per year. Q in and Q rej are decreasing 
since the total heat is decreasing as shown in Figure 5(b) for one ASC; while T h is maintained as constant 
as that at the BOM, T c decreases a little as shown in Figure 5(c), which results in a higher ASC efficiency. 
The piston amplitude at EOM is approximately 91.1 percent of that at the BOM. 

For the case of starting from the minimum fuel load at BOM, the ASRG DC power output is 121-We 
with the ASRG efficiency of 28.4 percent at EOM, compared with a 140-We DC power output with 
28.6 percent ASRG efficiency at BOM in deep space, which is referred as the current best estimation 
(CBE) ASRG power output in (Ref. 7). 

The corresponding ASRG power output with the second ASC control strategy( fixed piston amplitude) 
is presented in Figure 6, showing ASRG DC power output drops 30.5 We within the 17-year period 
starting from the minimum fuel load at BOM. Q in and T h are decreasing with YOM, and Q rej and T c 
remain almost constant. The ASC efficiency drop is more significant; thus power output drops more. At 
EOM, the ASRG DC power output is 109.5-We with the efficiency of 25.7 percent. 

For the case of starting from the maximum fuel load at BOM, the ASRG output power drop is very 
similar to that starting at the minimum fuel load. Also we can see that the DC power output is decreasing 
linearly with Y OM since the fuel decay is almost linear. 


YOM GPHS Fuel loading 
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Figure 5. — ASRG performance under YOM with fixing Th. 
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Figure 6. — ASRG performance under YOM with a fixed piston amplitude. 

3.2 ASRG Performance Characteristic Under Different Ambient Sink Temperature 

For different missions, the environment could be very different. Here, the ambient sink temperature 
T a varying from 4 to 300 K as shown in (Ref. 7) is investigated. The equivalent sink temperature of 200 K 
is used for the environment on Mars. For a Venus flyby, the sink temperature can be 300 K. It is assumed 
here that heat only radiates between the ASRG housing and the ambient. The GPHS fuel load is assumed 
to be 244, 250, and 256 W for minimum, nominal, and maximam cases at BOM, respectively. The two 
ASC control strategies are also investigated here. 

Figure 7 shows how the ASRG power output varies with the sink temperature for the first ASC 
control strategy. When the sink temperature gets higher — that is, a warmer environment — Q m increases a 
little because insulation loss decreases shown in Figure 7(b) for one ASC. If the 77, is fixed, 77 increases 
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while T a gets higher as shown in Figure 7(c), which results in a lower ASC efficiency. Thus the power 
output will drop. For the minimum fuel load case, at T a = 200 K, we can see a 5.2 We power drop from T a 
= 4 K, while at T a = 300 K, we can see that 21 -We drop and the ASRG efficiency becomes 24.4 percent. 
In order to maintain T h constant, the piston amplitude at T a = 300 K becomes 97.8 percent of that at 
T a = 4 K. 

Figure 8 shows the corresponding result using the second ASC control strategy. For the minimum 
fuel load case, the ASRG DC power output drops 7 We when T a = 200K, and 27 We drop when 
T a = 300 K. If the piston amplitude is fixed as a constant, T h drops and T c increases while T a gets higher as 
shown in Figure 8(c), which results in a larger ASC efficiency drop. At T a = 300 K, the ASRG efficiency 
becomes 23.2 percent. 

It shows that the first ASC control strategy can achieve higher ASC efficiency when the ambient gets 
warmer, thus the power drop becomes less than that when using the second ASC control strategy. 
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Figure 7. — ASRG performance under varying sink temperature T a with fixing Th. 
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Figure 8. — ASRG performance under varying sink temperature T a with a fixed piston amplitude. 


4.0 Conclusions 

A one-dimensional steady-state mathematical thermal power model of the ASRG has been presented. 
How the ASC represented by Sage performance map is coupled in the ASRG thermal circuit has been 
explained. The ASRG performance characteristics under different environments and GPHS fuel load 
decay with Y OM have been shown. Further validation of the one-dimensional model with a three- 
dimensional finite element ASRG thermal-power model for the flight unit or EU will be performed in 
future work. 
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Appendix — Symbols List 


A i housing surface area 

A 2 heat transfer area 

A p piston amplitude 

/ frequency of the Stirling engine 

h convection heat transfer coefficient 

m mass of the working fluid 

P the power output of Stirling engine 

p m mean pressure of the ASC 
Q in heat going to the ASC 

Qiossj heat loss from the GPHS through insulation material to the ambient 

Qioss 2 heat loss from the ASC hot end through insulation material to the ambient 

Qhssj heat loss through conduction inside the ASC 

Q rej heat rejected from the cold end of the ASC 

Qtotai total heat generated from the GPHS 

Qwork heat going to do the work 

R gas constant 

Ri contact and conduction resistances between GPHS and housing wall 

R 2 contact and conduction resistances between GPHS and heat collector of the ASC 

/C, contact and conduction resistances between the ASC hot end and housing wall 

R/[ conduction resistance between the cold end of the ASC and housing wall 

R' 4 convection and/or radiation resistances between housing wall and the ambient environment 

R s conduction resistance between the hot and cold ends of the ASC 

R () h resistance at the hot end of the ASC 
R 6c resistance at the cold end of the ASC 

R‘ 6 h resistance at the hot end of an idealized Stirling engine 

R'e c resistance at the cold end of an idealized Stirling engine 

i? 6 u h resistance at the hot end based on the West equation 

R” c resistance at the cold end based on the West equation 

T a equivalent ambient sink temperature 

T h temperature at the hot end of the ASC 

T c temperature at the cold end of the ASC 

T s temperature of the GPHS 

T w housing wall temperature 

V D displacer swept volume 

V p piston swept volume 

W the mechanical work from Stirling engine 

w„ West number 

e surface emissivity 

a Stefan-Boltzmann constant 
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